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Abstract 

We report the results of HCN(J=4-3) and HCO"'"( J=4-3) observations of two lu- 
minous infrared galaxies (LIRGs), NGC 4418 and Arp 220, made using the Atacama 
Submillimeter Telescope Experiment (ASTE). The ASTE wide-band correlator pro- 
vided simultaneous observations of HCN(4-3) and HCO^(4-3) lines, and a precise 
determination of their flux ratios. Both galaxies showed high HCN(4-3) to HC0"'"(4- 
3) flux ratios of >2, possibly due to AGN-related phenomena. The J = 4-3 to J = 
1-0 transition flux ratios for HCN (HCO"'') are similar to those expected for fully ther- 
malized (sub-thermally excited) gas in both sources, in spite of HCN's higher critical 
density. If we assume coUisional excitation and neglect an infrared radiative pumping 
process, our non-LTE analysis suggests that HCN traces gas with significantly higher 
density than HCO"*". In Arp 220, we separated the double-peaked HCN(4-3) emis- 
sion into the eastern and western nuclei, based on velocity information. We confirmed 
that the eastern nucleus showed a higher HCN(4-3) to HCN(l-O) flux ratio, and thus 
contained a larger amount of highly excited molecular gas than the western nucleus. 
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galaxies — galaxies: individual (NGC 4418 and Arp 220) 
1. Introduction 

Luminous infrared galaxies (LIRGs) radiate the bulk of their large luminosities (L ^ 
IO^^Lq) as infrared dust emission (Sanders & Mirabel 1996). Their large infrared luminosities 
indicate that (1) powerful energy sources are present, hidden behind dust; (2) energetic radiation 
from the hidden energy sources is absorbed by the surrounding dust; and (3) the heated dust 
grains re-emit this energy as infrared thermal radiation. To understand this LIRG population, 
it is essential to clarify the hidden energy sources, whether starbursts (i.e., nuclear fusion inside 
stars) are dominant, or active galactic nuclei (AGNs; active mass accretion onto a central 
compact supermassive black hole [SMBH] with >10^Mq) are also energetically important. 

Unlike AGNs surrounded by a torus-shaped (toroidally-shaped) dusty medium, which 
are classified optically as Seyferts (Veilleux & Osterbrock 1987), LIRGs have a large amount of 
concentrated molecular gas and dust in their nuclei (Sanders & Mirabel 1996). For this reason, 
the putative compact AGNs can easily be buried (i.e., obscured in virtually all directions), 
making the optical detection of AGN signatures very difficult. However, investigating the 
energetic importance of such optically elusive buried AGNs is crucial to understand the true 
nature of the LIRG population. 

A starburst (nuclear fusion) and a buried AGN (mass accretion onto a SMBH) have very 
different energy-generation mechanisms. Specifically, while UV is the predominant energetic 
radiation in a starburst, an AGN emits strong X-ray emission in addition to UV. The radiative 
energy generation efficiency of a nuclear fusion reaction in a normal starburst is only ~0.5% of 
Mc^ (M is the mass of material used in the nuclear fusion reaction). Thus, the emission surface 
brightness of a starburst region is modest and has both observational (Werner et al. 1976; Soifer 
et al. 2000) and theoretical (Thompson et al. 2005) upper limits (~10^^Lq kpc"^). An AGN, 
however, achieves high radiative energy generation efficiency (6-42% of Mc^, where M is the 
mass of accreting material; Bardeen 1970; Thorne 1974), producing a very high emission surface 
brightness with >10^^Lq kpc~^ (Soifer et al. 2000). These differences between a starburst and 
an AGN could create differences in the properties of the molecular gas and dust around the 
energy sources, allowing us to distinguish the hidden energy sources if, for example, different 
emission-line flux ratios are observed from the surrounding material. 

Molecular gas emission lines in the sub-millimeter and millimeter wavelength ranges are 
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a potentially effective tool for probing this distinction, because the effects of dust extinction 
are so small in these wavelength ranges that signatures of AGNs deeply buried in gas and 
dust are detectable. It was found observationally that HCN(J=l-0) emission (rest frequency 
I'rest = 88.632 GHz or rest wavelength Arcst = 3.385 mm) is systematically stronger, relative 
to HCO"^(J=1-0) (i^rcst = 89.189 GHz or Arest = 3.364 mm), in AGN-dominated galaxy nuclei 
than in starburst galaxies (Kohno 2005; Krips et al. 2008). The trend of a strong HCN(l-O) to 
HCO"*"(1-0) emission ratio was further confirmed in luminous buried AGN candidates (Imanishi 
et al. 2004; Imanishi et al. 2006b; Imanishi & Nakanishi 2006; Imanishi et al. 2007b; Imanishi 
et al. 2009). These results demonstrated that this HCN(1-0)/HCO^(1-0) method functions as 
a diagnostic tool to distinguish the hidden energy sources of dusty galaxies. Extension to more 
distant, faint LIRGs is expected to play an important role for the comprehensive understanding 
of the LIRG population in the forthcoming Atacama Large Millimeter /submillimeter Array 
(ALMA) era. However, during the first operation phase of ALMA, band 3 (84-116 GHz or 
2.59-3.57 mm) is the longest available wavelength coverage, so that this HCN(1-0)/HCO"^(1-0) 
method can be applied only to LIRGs a.t z< 0.05, hampering further systematic investigation of 
LIRGs. Energy diagnostic methods using other transition lines and/or other molecular species 
must be exploited. 

HCN( J=4-3) (i/rest = 354.505 GHz) and HCO+(4-3) (u.^st = 356.734 GHz) lines exist in 
Earth's 350 GHz atmospheric window, and can be one of the main frequency ranges exploited 
with ALMA. Because the frequency coverage of ALMA band 7 is 275-373 GHz, HCN(4-3) and 
HCO"*"(4-3) lines can be studied in LIRGs aX z > 0.05. Thus it is important to investigate 
whether AGNs show high HCN(4-3) to HCO+(4-3) flux ratios, similar to the high HCN(l-O) 
to HCO+(1-0) flux ratios. 

In this paper, we present the results of simultaneous HCN(4-3) and HCO"''(4-3) obser- 
vations of two nearby, well-studied LIRGs, NGC 4418 and Arp 220, both of which show high 
HCN(l-O) to HCO+(1-0) flux ratios (Imanishi et al. 2004; Imanishi et al. 2007b). Throughout 
this paper, we adopt Hq = 75 km s~^ Mpc~^, Qm = 0.3, and Q\ = 0.7, to be consistent with 
our previously published papers. 

2. Targets 

NGC 4418, at z = 0.007, is a nearby LIRG with an infrared luminosity of Lir ~ 9 
X 1O^°L0 (Table 1). Although it shows no clear Seyfert signature in the optical spectrum 
(Armus, Heckman, & Miley 1989; Lehnert & Heckman 1995), it is considered one of the closest 
and strongest buried AGN candidates, based on the AGN-like infrared spectral shape (Roche 
et al. 1991; Dudley & Wynn- Williams 1997; Spoon et al. 2001) and the high HCN(l-O) to 
HCO"*"(1-0) flux ratio in the millimeter wavelength range (Imanishi et al. 2004). 

Arp 220 (z = 0.018) is one of the best-studied nearby ultraluminous infrared galaxies 
(LiR > IO^^Lq; Sanders et al. 1988). It has two nuclei, the eastern (Arp 220 E) and western 
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(Arp 220 W) nuclei, separated by ~1" (Scoville et al. 2000; Soifer et al. 2000). The combined 
optical spectrum of both nuclei is classified as a LINER (i.e., non-Seyfert) (Veilleux et al. 1999). 
Starburst activity has been detected, but was energetically insufficient to account quantitatively 
for the observed infrared luminosity (Imanishi et al. 2006a; Armus et al. 2007; Imanishi et al. 
2007a), indicating an energy source deeply buried in Arp 220's nuclear core (Spoon et al. 
2004; Gonzalez- Alfonso et al. 2004). The energy source in Arp 220 W has a very high emission 
surface brightness, for which a buried AGN is a plausible (if not definite) explanation (Downes 
& Eckart 2007; Sakamoto et al. 2008). 

3. Observations and Data Reduction 

The HCN(4-3) and HCO+(4-3) observations of NGC 4418 and Arp 220 were made using 
the Wideband and High dispersion Spectrometer system with FFX correlator (WHSF) (Okuda 
& Iguchi 2008; Iguchi & Okuda 2008) on the Atacama Submillimeter Telescope Experiment 
(ASTE), a 10-m antenna located at Pampa La Bola in the Atacama Desert of Chile, at an 
altitude of 4800 m (Ezawa et al. 2004; Ezawa et al. 2008). Because the WHSF has 4-GHz 
frequency coverage in each side band, both the HCN(4-3) (z-'rcst = 354.505 GHz) and HC0+(4- 
3) (t'rest = 356.734 GHz) lines were simultaneously observable. Figure 1 shows the raw ASTE 
WHSF spectra of NGC 4418 and Arp 220. 

The observations were made in 2008 May, remotely from an ASTE operation room 
in the Mitaka campus of the National Astronomical Observatory of Japan, using the network 
observation system N-C0SM0S3 developed by Kamazaki et al. (2005). Table 2 summarizes the 
detailed observation log. The weather conditions were excellent. Typical system temperatures 
were 180-300 K in the single side band. We employed the position-switching mode, so that the 
target position and off-sky position, 4 arcmin away from the target, were switched after each 
10 s exposure. Telescope pointing was checked every 2 h. The total net on-source integration 
times were 200 min and 155 min for NGC 4418 and Arp 220, respectively. Flux calibration 
was made every observation day. By changing the frequency setting of the WHSF, we observed 
a strong emission line in W28 (IRC10216) at the frequency of the red-shifted HCN(4-3) and 
HCO^(4-3) lines of NGC 4418 (Arp 220), to estimate the possible sensitivity variation along 
the WHSF bandpass. This variation was found to be <20%. The half-power beam width of 
the ASTE 10-m dish was 22" at ~ 345 GHz. At this frequency, the main beam efficiency is 
0.6, and 1 K in antenna temperature corresponds to 78.5 Jy. 

Data reduction was carried out with NEWSTAR, a package developed at Nobeyama 
Radio Observatory. A fraction of poor-quality data was flagged. To estimate the strength of 
HCN(4-3) and HCO+(4-3) emission, we subtracted baselines from the obtained raw spectra. 
We first divided the spectra into two parts at the center, the lower-frequency HCN(4-3) part 
and higher-frequency HCO^(4-3) part. The baseline was estimated independently in each part, 
based on the data points free from emission lines. 

4 



For NGC 4418, ten data sets with 20 min net on-source exposure time for each were 
obtained. The behavior of the basehne was extremely good, and the HCN(4-3) emission hne was 
easily recognizable in individual data sets. First-order linear baselines were adopted both for 
the HCN(4-3) and HCO^(4-3) emission lines, and the baseline-subtracted individual spectra 
were summed to obtain final spectra. 

For Arp 220, final data consist of seven data sets with 20 min net on-source exposure 
and one data set with 15 min net on-source exposure. For the HCN(4-3) line, first-order linear 
baselines were adopted in individual data sets, and the baseline-subtracted individual spectra 
were summed. The HCN(4-3) emission line was detected in the final spectrum (see Appendix 
for more detail). Although we tried several baseline fits, the HCN(4-3) emission line was 
always visible, and its profile was virtually unchanged in the final spectrum, under reasonable 
baseline choices. For the HCO^(4-3) line, when we applied a first-order linear baseline, winding 
patterns remained in the final spectrum, and HCO"''(4-3) line was not clearly detected. We 
thus attempted a second-order baseline, but the HCO"'"(4-3) emission line remained undetected 
in the final spectrum. 

4. Results 

Figure 2 presents the final spectra around the HCN(4-3) and HCO"'"(4-3) emission lines, 
after baseline subtraction. In both NGC 4418 and Arp 220, HCN(4-3) emission was clearly 
detected, with peak antenna temperatures of 7-15 mK, but HCO"'"(4-3) emission was not. For 
NGC 4418, HCN(4-3) and HCO"'"(4-3) emission line spectra are published here for the first 
time. 

For NGC 4418, we fit the detected HCN(4-3) emission with a single Gaussian profile. For 
Arp 220 HCN(4-3) emission, a double-peaked profile was clearly seen. Thus we fit it with a two- 
component Gaussian profile. Table 3 summarizes the fitting results. The estimated HCN(4-3) 
fluxes were 220 Jy km s~^ and 410 Jy km s~^ for NGC 4418 and Arp 220, respectively, with 
the uncertainties of ~20%. For the undetected HCO"*"(4-3) emission, we adopted the lowest 
plausible continuum level in the spectra, and estimated conservative upper limits, by counting 
the possible flux excess above the adopted continuum levels. Table 3 summarizes these upper 
limits for the HCO"''(4-3) emission. 

For Arp 220, HCN(4-3) measurements have been made previously by two groups 
(Wiedner et al. 2002; Greve et al. 2009). Our ASTE measurement of the HCN(4-3) flux 
(410±82 Jy km s~^) fell between the previously reported values of 260 Jy km s~^ (Wiedner 
et al. 2002) and 587±118 Jy km s~^ (Greve et al. 2009). The velocity difference in the two 
components in the double-peaked HCN(4-3) profile was ~300 km s~^ (Table 3), similar to 
that of HCN(l-O) emission (Imanishi et al. 2007b), further strengthening the detection of the 
HCN(4-3) emission line in Arp 220. Based on the argument of Imanishi et al. (2007b), we 
ascribe the lower-velocity (blue-shifted) and higher-velocity (red-shifted) components to the 
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Arp 220 W and E nuclei, respectively. 

The HCN(4-3) to HCO+(4-3) flux ratios ^ were >2.7 and >3.0 for NGC 4418 and Arp 
220, respectively, where the ratio for Arp 220 used the combined emission from both nuclei. 
The non-detection of HCO+(4-3) emission lines hampered the spectral separation into Arp 220 
E and W and rendered independent estimates of the line flux ratio for these nuclear components 
difficult. HCN(4-3)/HCO"'"(4-3) flux ratios significantly higher than unity were confirmed in 
both galaxies. 

We emphasize that our simultaneous HCN(4-3) and HCO"'"(4-3) observations us- 
ing ASTE WHSF are crucial for obtaining rehable HCN(4-3) to HCO+(4-3) flux ratios. 
Submillimeter and millimeter observations still present non-negligible uncertainties in abso- 
lute flux calibrations. In fact, even for measurements obtained in the last several years, the 
HCN(4-3) and HON (1-0) fluxes of Arp 220 differed by a factor of 2 at maximum as observed 
by different groups with different telescopes under different weather conditions and using dif- 
ferent calibration methods (Wiedner et al. 2002; Imanishi et al. 2007b; Gracia-Carpio et al. 
2008; Greve et al. 2009). This level of absolute flux calibration uncertainty can be fatal to 
an investigation of the HCN(4-3) to HCO"'"(4-3) flux ratios because we are discussing the 
difference of the ratios by about a factor of 2. By simultaneously observing HCN(4-3) and 
HCO^(4-3) lines, the possible absolute flux calibration uncertainty does not propagate to the 
HCN(4-3)/HCO"''(4-3) flux ratios, which are then dominated only by statistical noises and 
fltting errors. 

5. Discussion 

5. 1 . HCN(4 -3) to HCO+ (4 -3 ) flux ratios 

In this subsection, we explore HCN(4-3) and HCO+(4-3) emission from NGC 4418 and 
Arp 220, in relation to HCN(l-O) and HCO+(1-0). For this purpose, we used our ASTE 
WHSF and Nobeyama Millimeter Array (NMA) data (Imanishi et al. 2004; Imanishi et al. 
2007b), because (1) simultaneous HCN and HCO"*" observations by both telescopes make the 
comparison of HCN to HCO+ flux ratios rehable, and (2) HCN(l-O) and HCO+(1-0) absolute 
fluxes in the NMA data have been conflrmed to agree with measurements from other single-dish 
telescopes (Imanishi et al. 2004; Imanishi et al. 2007b). The NMA measurements of HCN(l-O) 
and HCO"'~(1-0) fluxes are summarized in Table 4. 

For NGC 4418, we obtained an HCN(4-3) to HCO+(4-3) flux ratio > 2.7, conflrming 
the high HCN to HCO+ flux ratio previously found for the J = 1-0 transition line (~1.8; 
Imanishi et al. 2004). The HCN(4-3) to HCN(l-O) flux ratio was estimated to be 20. This is 

^ By definition, flux is proportional to brightness-temperature x v^, where v is frequency. Because HCN (4-3) 
(frest = 354.505 GHz) and HCO+(4-3) (t'rest = 356.734 GHz) have very similar frequencies, their flux ratios 
are virtually identical to their brightness-temperature ratios. In this paper, we use a flux ratio. 
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comparable to the scaling (~16) expected for fully thermalized optically thick gas (Riechers 
et al. 2006). The HCO+(4-3) to HCO+(1-0) flux ratio was <15. 

For Arp 220, the HCN(4-3) and HCO+(4-3) flux ratio was estimated to be >3.0. The 
J = 4-3 to J = 1-0 flux ratio for Arp 220 total emission was 12 and <8 for HCN and HCO+, 
respectively. When we separated the emission from individual nuclei, the HCN(4-3) to HCN(1- 
0) flux ratios were 21 and 9 for Arp 220 E and W nuclei, respectively. The ratio is higher in 
Arp 220 E than Arp 220 W, further supporting the previous argument that a larger fraction of 
highly excited molecular gas is present in Arp 220 E (Greve et al. 2009). 

The critical density of HCN is 2.3 x 10^ cm~^ for J = 1-0 and 8.5 x 10^ cm"^ for J = 
4-3. That of HCO+ is 3.4 x 10"^ cm'^ for J = 1-0, and 1.8 x 10*^ cm'^ for J = 4-3 (Table 3 
of Greve et al. 2009). Because the critical density of HCO"'" is a factor of ~5-7 smaller than 
HCN for both the J = 1-0 and J = 4-3 transitions, HCO"'" can be coUisionally excited more 
easily than HCN in single-phase molecular gas. Thus, the stronger HCN emission in NGC 4418 
and Arp 220 requires (1) an HCN abundance that is significantly higher than HCO"*", and/or 
(2) the presence of some physical mechanism that enhances HCN emission. 

Regarding the first HCN abundance enhancement scenario, Yamada et al. (2007) con- 
cluded, based on three-dimensional line transfer simulations, that HCN abundance must be an 
order of magnitude higher than HCO"*" to account for the observed high HCN to HCO^ flux 
ratio (~2). While the typical HCN to HCO"^ abundance ratio in our Galactic molecular cloud 
is ~1 (Blake et al. 1987; Pratap et al. 1997; Dickens et al. 2000), HCN overabundance, relative 
to HCO"*", by an order of magnitude, is predicted in some parameter range for dense molecular 
gas illuminated by an X-ray emitting energy source (Meijerink & Spaans 2005). Since an AGN 
is a much stronger X-ray emitter than starburt activity, this HCN overabundance scenario 
is a possibility for the observed high HCN(4-3)/HCO"'"(4-3) flux ratios in the buried AGN 
candidates, NGC 4418 and Arp 220. 

Strong HCN emission in luminous buried AGN candidates may be explained by the 
infrared radiative pumping scenario (Aalto et al. 1995; Garcia-Burillo et al. 2006; Guelin et 
al. 2007; Weiss et al. 2007; Aalto et al. 2007a). In an AGN, the emission surface brightness is 
much higher than in a starburst, so that the surrounding dust is heated to a high temperature 
(several hundred K) and produces strong mid-infrared 10-20 //m continuum emission. The HCN 
molecule has a line at infrared 14.0 /im. Molecular gas around an AGN can be vibrationally 
excited by absorbing these infrared 14.0 ^m photons, and through the subsequent cascade 
process, can enhance HCN rotational lines in the sub-millimeter and millimeter range. 

Guelin et al. (2007) suggested that this infrared pumping scenario may work also for 
HCO+ and HNC, because they have lines at 12.1 yum (HCO"^) and 21.7 yum (HNC). However, 
for this infrared pumping scenario to affect particular molecules, infrared photons at the wave- 
lengths of these molecular lines must be absorbed. Namely, absorption features must be de- 
tected in infrared spectra at 14.0 /xm (HCN), 12.1 /xm (HCO+), and 21.7 /im (HNC). While 
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the infrared HCN 14.0 /im absorption feature is clearly detected in highly obscured ultralumi- 
nous infrared galaxies, including NGC 4418 and Arp 220 (Lahuis et al. 2007; Veilleux et al. 
2009), absorption features at HCO+ 12.1 fim and HNC 21.7 /im are not (Veilleux et al. 2009). 
Therefore, we see observational signatures that this infrared radiative pumping scenario could 
be at work for HCN, but not for HCO+ or HNC. This may explain the strong HCN emission 
found in luminous AGN candidates. 

Since both the HCN overabundance and infrared radiative pumping scenarios could work 
more effectively in an AGN than starburst activity, the enhanced HCN(4-3) emission may 
originate in putative buried AGNs in NGC 4418 and Arp 220. Three-dimensional modelings 
which realistically incorporate actual physical parameters of molecular gas around an AGN are 
needed for more detailed quantitative comparison with observations. 

5.2. Molecular gas properties 

For NGC 4418, the HCN(4-3) to HCN(l-O) flux ratio (20) was higher than the HC0+(4- 
3) to HCO^(l-O) flux ratio (<15). For Arp 220 as well, the former ratio (12) was higher than 
the latter ratio (<8). Thus, our ASTE observations suggest that J = 4-3 transition lines are 
excited more efficiently in HCN than HCO"'", in spite of the higher critical density of HCN(4- 
3) compared to HCO"^(4-3) (a factor of ~5 [§5.1]). As mentioned in §5.1, in a single-phase 
molecular gas, molecular species with lower critical density are more easily coUisionally excited. 
Thus, HCO"'" should be excited more than HCN. Some mechanisms are needed to explain the 
strong HCN(4-3) emission. Examples include a multiple phase molecular gas (Greve et al. 
2009), and the infrared radiative pumping scenario. The second scenario will be investigated 
elsewhere, using a three-dimensional model (Yamada et al. 2009, in prep). Here, we explore 
the first possibility. Namely, we consider only collisional excitation and neglect the infrared 
radiative pumping process. 

We used RADEX (van der Tak et al. 2007) to investigate molecular gas properties that 
could account for our ASTE and NMA results. The RADEX software adopts the non-LTE 
analysis of molecular line spectra, solving radiative transfer in simple geometry gas. We varied 
the H2 number density and kinetic temperature from lO'^-lO^ cm~'^ and 10-200 K, respectively. 
A uniform sphere with a line width of 200 km s~^ was assumed, based on the HCN(4-3) emission 
measurements (Table 3). X-ray observations suggested that both NGC 4418 and Arp 220 have 
highly-obscured, Compton thick (Nh > lO^'^ cm~^) AGNs (Maiolino et al. 2003; Iwasawa et 
al. 2005). Taking HCN and HCO+ abundances relative to H2 of 10^^-10"^ (Blake et al. 
1987; Dickens et al. 2000), column densities with 10^^ cm~^ were adopted for both HCN and 
HCO"^. The T = 2.73 K background emission was also added for the calculation. Figure 3 
shows the HCN(4-3) to HCN(l-O) and HCO+(4-3) to HCO+(1-0) fiux ratios in Jy km s-\ 
as a function of molecular gas number density and kinetic temperature. In the 10-200 K 
temperature range, the molecular gas number density was estimated to be uh >10^ cm""^ to fit 
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the observed HCN(4-3) to HCN(l-O) flux ratios for NGC 4418 (~20) and Arp 220 (~12). 

Now, the HCO+(4-3) to HCO+(1-0) Hux ratios were <15 and <8 for NGC 4418 and 
Arp 220, respectively. For T > 25 K (>10^'^ K), molecular gas with a number density uh 
<10^ cm~^ reproduced the observed HCO+ flux ratios. Thus, assuming collisional excitation, 
we confirmed that in both NGC 4418 and Arp 220, HCN emission traces gas with a number 
density significantly higher than that traced by HCO"'" ^, as previously argued for Arp 220 
based on the measurements of HCN(4-3), HCO"*~(4-3), and other transition lines (Greve et al. 
2009). Even given HCN and HCO"'" column densities of 10^^ cm~^ and a line width of 300 km 
s~^, our main conclusion is unchanged. 

Three-dimensional hydrodynamic simulations of galactic disks have shown that the frac- 
tion of molecular gas with uh > 10^ cm~'^ is miniscule (Wada & Norman 2007). Yet, the strong 
HCN (4-3) emission line detected in NGC 4418 and Arp 220 suggests the presence of a large 
amount of such high-density (uh > 10^ cm~^) molecular gas, if the HCN(4-3) line is coUisionally 
excited. It is likely that the high nuclear concentration of molecular gas in LlRGs increases the 
fraction of high-density molecular gas. 

6. Summary 

We presented the results of simultaneous HCN (4-3) and HCO"*" (4-3) observations of two 
luminous infrared galaxies, NGC 4418 and Arp 220, using ASTE WHSF. Our main conclusions 
are as follows: 

1. Strong HCN(4-3) emission was detected in both NGC 4418 and Arp 220, but HCO+(4-3) 
emission provided only upper limits. 

2. In both sources, we found HCN(4-3) to HCO"*" (4-3) flux ratios of >2, even higher than 
the HCN(l-O) to HCO"'" (1-0) emission hue flux ratios. The high HCN/HCO+ flux ratios 
could be explained by HCN overabundance and/or infrared radiative pumping scenarios, 
both of which may work more effectively in an AGN than starburst activity. 

3. HCN(4-3) to HCN(l-O) flux ratios were higher than HCO+(4-3) to HCO+(1-0) flux ratios 
in both NGC 4418 and Arp 220, suggesting that higher J-transitions were excited more 
efficiently in HCN than HCO"*". If we neglect an infrared radiative pumping mechanism, 
and consider only collisional excitation, these results are incompatible with a one-zone 
molecular gas model, because the critical density of HCN is higher than HCO'^ at these 
transitions. A multi-phase molecular gas in which HCN selectively probes higher-density 
molecular gas than HCO"*" is required. We applied RADEX models, and found that in 
both NGC 4418 and Arp 220, gas with a number density uh > 10^ cm~^ and uh < 10^ 

2 For NGC 4418, the small HCN(3-2) to HCO+(3-2) flux ratio (less than unity) reported by Aalto et al. 
(2007b) is not compatible with our model. The HCN(3-2) flux itself is also significantly lower than that 
expected for thermalized, optically thick gas. We have no clear interpretation, as long as only collisional 
excitation is considered. 
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cm ^ reproduced the observed HCN and HCO"*" J = 4-3 to 1-0 flux ratios, respectively, 
in the 25-200 K temperature range. 
4. For Arp 220, HCN(4-3) to HCN(l-O) flux ratio is higher in Arp 220 E than in Arp 220 
W, suggesting that Arp 220 E contains a larger amount of highly excited molecular gas. 

We thank T. Okuda for the support of our ASTE WHSF observing run, and the anony- 
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in collaboration with University of Chile, and Japanese institutes including University of 
Tokyo, Nagoya University, Osaka Prefecture University, Ibaraki University, and Hokkaido 
University. Observations with ASTE were in part carried out remotely from Japan by us- 
ing NTT's GEMnet2 and its partner R&E (Research and Education) networks, which are 
based on AccessNova collaboration of University of Chile, NTT Laboratories, and NAOJ. 
M.I. is supported by Grants-in-Aid for Scientific Research (19740109). This study utilized 
the NASA/IPAC Extragalactic Database (NED) operated by the Jet Propulsion Laboratory, 
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Appendix 1. Analysis of the Arp 220 HCN(4— 3) emission line 

We explain in detail how the HCN(4-3) emission line of Arp 220 becomes detectable in 
our analysis. We took eight independent data sets. Single-order linear baselines were employed 
in individual data sets (Figure 4) and eight baseline-subtracted spectra were summed to obtain 
a final HCN(4-3) spectrum. The HCN(4-3) emission was clearly detected in the final spectrum 
of Arp 220, as seen in Figure 2. 
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Table 1. Detailed information on the observed LIRGs 







J 12 


J25 


f 

J 60 


f 

J 100 


Inc T,TT~, lino' T,TT~> / 7^,^ 1 






(Jy) 


(Jy) 


(Jy) 


(Jy) 


(ergs s~-^) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


NGC 4418 


0.007 " 


0.9 


9.3 


40.7 


32.8 


44.5 (10.9) 


Arp 220 


0.018 


0.48 


7.92 


103.33 


112.40 


45.7 (12.1) 



Notes. 

Col.(l): Object name. 
Col. (2): Redshift. 

Cols.(3)-(6): fi2, f25, feo; and fioo are IRAS FSC fluxes at 12, 25, 60, and 100 /im, 
respectively. 

Col. (7): Decimal logarithm of the infrared (8—1000 /im) luminosity in ergs s~^ calculated 
as follows: Ljr = 2.1 x lO^'^x D(Mpc)2 x (13.48 x /12 + 5.16 x /as + 2.58 x /so + /loo) ergs 
s~^ (Sanders & Mirabel 1996). The values in parentheses are the decimal logarithms of the 
infrared luminosities in units of solar luminosities. 
1 arcsec corresponds to 130 pc. 

^: 1 arcsec corresponds to 340 pc. 



Table 2. Observing log of ASTE observations 



Object 


Observing Date 


Central 


Pointing 


Flux 




(UT) 


frequency 


calibrator 


calibrator 


(1) 


(2) 


(3) 


(4) 


(5) 


NGC 4418 


2008 May 22, 23, 24, 27 


353.15 


V Hya, IRC10216 


W28 


Arp 220 


2008 May 24, 27 


349.21 


RX Boo 


IRC10216 



Notes. 

Col.(l): Object name. 

Col. (2): Observation date in UT. 

Col. (3): Central frequency of ASTE WHSF used for the observations. 
Col. (4): Pointing calibrator object name. 
Col. (5): Flux calibrator object name. 
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Table 3. Properties of HCN(4-3) and HCO+(4-3) emission lines 



Object 


Line 


LSR velocity 


Line width 


Flux 


Flux 


Luminosity 


Ratio 






[km s J 


[km s J 


[K km s^^] 


[Jy km s"^] 


lU [K km s pc J 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


NGC 4418 


HCN(4-3) 


2060 


335 


4.7 


220 


4.7 


>2.7 




HCO+(4-3) 






<1.7 


<80 


<1.7 




Arp 220 


HCN(4-3) 


5265, 5560 


220,200 


8.7 


410 


57 


>3.0 










3.4 (W) + 5.3 (E) 


160 (W) + 250 (E) 


22 (W) + 35 (E) 






HCO+(4-3) 






<2.9 


<135 


<19 





Notes. 

Col.(l): Object name. 

Col. (2): HCN(4-3) or HCO+(4-3) line. 

Col.(3): LSR velocity {vopt = {f - 1) x c} of the HCN(4-3) and HCO+(4-3) in [km 
Col. (4): Line width in FWHM of the HCN(4-3) and HCO+(4-3) emission lines in [km 

s-i]. 

Col. (5): Flux of the HCN(4-3) and HCO^(4-3) emission lines in main beam temperature 
in [K km s~^]. For ASTE, the main beam efficiency is ~ 0.6 at ~ 345 GHz. For the HCN(4- 
3) line of Arp 220, individual contributions from Arp 220 E and W are also shown. The flux 
uncertainty is estimated to be <20%, based on the combination of the empirical one of ASTE 
WHSF (10-15%) and systematic one of flux calibrators (~10%) (Wang et al. 1994). 

Col. (6): Flux of the HCN(4-3) and HCO+(4-3) emission lines in [Jy km s'^]. For ASTE, 
1 [K] in main beam temperature corresponds to 47 [Jy] at ~ 345 GHz. 

Col.(7): Luminosity of the HCN(4-3) and HCO+(4-3) emission lines in xlO^[K km s"^ 

pc^]. 

Col.(8): HCN(4-3)/HCO+(4-3) flux ratio (flux oc u"^ x brightness-temperature). 
Because both the HCN(4-3) and HCO^(4-3) data are taken simultaneously, the ratio is not 
affected by possible absolute flux calibration uncertainties in the ASTE observations (see §4). 
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Table 4. HCN and HCO+ data at J = 4-3 and J = 1-0 transition lines for NGC 4418 and Arp 220 



Object 


Line 


Flux 


Reference 


(1) 


(2) 


(3) 


(4) 


NGC 4418 


HCN(4-3) 


220±44 


This work 




HCO+(4-3) 


<80 


This work 




HCN(l-O) 


10.3 


Imanishi et al. (2004) 




HCO+(1-0) 


5.5 


Imanishi et al. (2004) 


Arp 220 


HCN(4-3) 


410±82 


This work 






250 (E), 160 (W) 


This work 






260 


Wiedner et al. (2002) 






587±118 


Greve et al. (2009) 




HCO+(4-3) 


<135 


This work 






106±23 


Greve et al. (2009) 






80 


Sakamoto et al. (2009) 




HCN(1~0) 


34 


Imanishi et al. (2007b) 






12 (E), 18 (W) 


Imanishi et al. (2007b) 




HCO+(1-0) 


18 


Imanishi et al. (2007b) 






5.5 (E), 6.5 (W) 


Imanishi et al. (2007b) 



Notes. 

Col.(l): Object name. 

Col. (2): HCN or HCO+ line. 

Col. (3): Flux in [Jy km s"^]. For Arp 220, individual fluxes from Arp 220 E and W are 
also shown. For HCN(l-O) and HCO+(1-0) fluxes from Arp 220 E and W, only the spatially 
unresolved core emission is extracted. 

Col. (4): Reference. 



16 



in 
o 



o 
6 

ro 
I 

O 
X 

in 



o - 
I 

O 

in 
I 



NGC 4418 (ASTE Raw Spectrum) 



HCN(4-3) 



HC0+(4-3) 



352 353 354 355 
observed frequency (GHz) 



in 
o 



. Arp '2I26 (AStE' Raw' Spectrum) 



o - 



m 
o 

6 
I 



HCN(4-3) 



HC0'^(4-3) 

i 




348 349 350 
observed frequency (GHz) 



351 



Fig. 1. ASTE WHSF wide-band (4 GHz) spectra of NGC 4418 and Arp 220, before baseline subtraction. 
The abscissa is observed frequency in [GHz] and the ordinate is antenna temperature in [K] . The spectrum 
is smoothed to a velocity resolution of 60 km s^^ and 50 km s^^ for NGC 4418 and Arp 220, respectively. 
The expected frequency of the HCN(4-3) and HCO+(4-3) lines arc indicated as arrows. HCN(4-3) and 
HCO"''(4-"3) emission lines arc simultaneously covered. 
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Fig. 2. HCN(4-3) and HCO+(4-3) spectra of NGC 4418 and Arp 220, after baseline subtraction. 
First-order linear fits were basically adopted as the baselines, except for the Arp 220 HCO"'"(4-3) line, for 
which a second-order fit was employed. The abscissa is velocity in [km s~^] and the ordinate is antenna 
temperature in [K]. Spectra of NGC 4418 and Arp 220 are smoothed to velocity resolutions of 60 km s~^ 
and 50 km s~^, respectively. Gaussian fits arc over-plotted as dashed lines for the detected HCN(4-3) lines 
in NGC 4418 and Arp 220. For the Arp 220 HCN(4-3) line, we ascribe the blue (lower velocity) and red 
(higher velocity) components to Arp 220 W (denoted as "W") and Arp 220 E ("E"), respectively, based 
on Nobcyama Millimeter Array interferometric HCN(l-O) and HCO+(1-0) data (Imanishi et al. 2007b) 
(see §4). The horizontal solid line indicates the zero flux level. 
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log (Temperature) log (Temperature) 



Fig. 3. HCN(4-3) to HCN(l-O) (Left) and HCO+(4-3) to HCO+(1-0) {Right) flux ratios in [Jy km s"!] 
as a function of gas number density and kinetic temperature, calculated with RADEX. The abscissa and 
ordinate are the decimal logarithms of kinetic temperature in [K] and H2 gas number density in [cm~"^], 
respectively. The HCN and HCO+ column density is assumed to be 10^^ cm~^ (see §5.2). A velocity 
width of 200 km s'^ was adopted (see §5.2). Contours are 0.1, 1, 5, 12, 20, 50, and 100 for HCN, and 0.1, 
1, 5, 8, 15, 50, and 100 for HCO+. The measured ratios for NGC 4418 and Arp 220 are indicated as thick 
lines. In the calculated parameter range, the line optical depth is 0-100 (i.e., reliable range in RADEX 
calculations) . 
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Fig. 4. Adopted linear baseline fits for individual unbinned HCN(4-3) spectra of Arp 220 (dashed lines), 
and smoothed spectra with velocity resolutions of 50 km s~^ (solid lines). Signatures of the HCN(4-3) 
emission lines are seen in most of the individual spectra. 
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